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ABSTRACT 

The  objective  of  this  study  was  to  develop  a  technique  of  predicting 
shock  shapes  and  wall  shock  traces  encountered  in  .secondary  injection 
into  a  conical  rocket  nozzle.  The  analysis  is  based  on  an  analogy  with 
Sakurai's  second  order  blast  wave  theory,  with  allowance  for  axial  vari¬ 
ation  in  the  free  stream  Mach  number  in  the  primary  nozzle.  Model  tests 
indicate  spreading  of  the  shock  adjacent  to  the  wall  slightly  higher 
than  expected.  A  procedure  for  shock  prediction  is  presented  based  on 
one  empirically  determined  constant.  Independence  of  side  forces  in 
perpendicular  secondary  injection  planes  under  conditions  of  shock 
intersection  was  experimentally  verified. 
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Local  radius  of  primary  nozzle 
Frontal  area  of  equivalent  blunt  body 
Constants 
Sonic  velocity 

P  - 

Pressure  coefficient,  ^ 

o  P  U 

Z  CO  00 

Pressure  coefficient  at  nose  of  equivalent  blunt  body 

Drag  coefficient  of  equivalent  blunt  body 

Exit  diameter  of  primary  nozzle 

Throat  diameter  of  primary  nozzle 

Equivalent  throat  diameter  of  Injected  flow 

Drag  force  on  equivalent  blunt  body 

Energy  release  per  unit  area  of  blast  surface 

Correction  factor  to  Z  accounting  for  variation  of 
Mach  number  along  the  primary  nozzle 

Constant  in  Sakural's  second  order  blast  wave  solution 

Primary  nozzle  radius  at  the  location  of  the  shock  apex 

Distance  from  primary  nozzle  exit  to  center  cf  injection 
port,  measured  along  the  nozzle  wall  in  a  diametral  plane 

Mass  flow  rate  of  Injectant 

Mass  flow  rate  of  primary  flow 

Mach  ntzmber 

Mach  number  of  primary  stream  at  shock  apex 
Pressure 
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Pressure  at  nose  of  the  olunt  body 
r.adius  of  equivalent  blunt  body 


u'ast  wave  radius 
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Radius  of  blast  wave  Inside  the  primary  nozzle 
Terta  in  Sakurai's  second  ordj,r  blast  wave  solution 


Velocity 


E  \(Ct+  1) 


it) 


Distance  mec‘ured  in  a  diametral  plane  along  the  primary 
nozzle  wall  from  the  apex  of  the  separation  line 

Distance  meai^^ired  in  a  diametral  plane  along  the  primary 
nozzle  wall  from  the  apex  of  the  shock  line 

Distance  downstream  of  the  shock  apex  along  a  line  passing 
through  the  shock  apex  and  parallel  to  the  primary  nozzle 
centerline 

Integer  in  Sakurai '  s  second  order  blast  wave  solution  o^O, 
1,2  corresponding  to  plane,  cylindrical,  and  spherical 
shocks,  respectively 

Angle  between  the  local  tangent  to  the  surface  of  the 
equivalent  blunt  body  and  the  undisturbed  primary  flow 

Value  of  p  at  the  nose  of  the  equivalent  blunt  body 

Axial  distance  measured  along  the  nozzle  centerline  from 
the  throat 

Shock  standoff  distance  upstream  of  the  nose  f  the  equiva¬ 
lent  blunt  body  measured  along  the  primary  nozzle  wall 

Ratio  of  heat  capacities 

Constants  in  Sakurai's  second  order  blast  wave  solution 

Azimuthal  angle  to  shock  trace  on  wall  referenced  to  a 
diametral  plane  passing  through  the  center  of  the  injection 
port.  The  angle  is  measured  clockwise  looking  forward. 

Density 


1 .  Introduction 

The  first  theoretical  treatment  of  secondary  injection  was  made 
by  Ferrari  [1]  in  a  study  of  the  forces  produced  by  a  lateral  jet  issuing 
from  a  cylindrical  body  into  a  supersonic  flow  past  the  body.  Ferrari 
used  Newtonian  theory  to  predict  the  trajectory  of  the  jet,  shock  stand¬ 
off  distances,  and  pressure  distributions;  shock  shapes  were  not  deter¬ 
mined  experimentally.  Amick  and  Hays  [2]  were  the  first  to  postulate  a 
secondary  injection  flow  model  on  the  basis  of  experimental  results.  In 
tests  of  secondary  injection  from  flat  plates  and  cylindrical  afterbodies, 
Amick  and  Hays  identified  separation  zones  and  interaction  shocks  by 
means  of  a  simple  flow  visualization  technique. 

The  first  attempt  to  correlate  secondary  injection  shocks  with  blast 
waves  was  made  by  Broadwell  [3]  who  utilized  first  order  blast  wave  theory 
to  predict  shock  shapes.  Broadwell 's  analysis  provided  a  good  framework 
for  data  correlation,  but  significantly  underpredicted  measured  side 
forces.  The  concept  of  treating  the  secondary  injection  jet  as  an  equiva¬ 
lent  blunt  body  was  first  put  forth  by  Zukoski  and  Spaid  [4]  in  a  study 
of  flat  plate  secondary  injection;  Schlieren  techniques  were  utilized  to 
determine  bow  shock  profiles  but  no  correlation  of  wall  shock  intersections 
were  made. 

A  significant  contribution  toward  formulation  of  a  flow  model  appli¬ 
cable  to  secondary  injection  in  a  rocket  nozzle  was  made  by  Charwat  and 
Allegre  [5] ;  these  experiments  identified  separation  and  shock  lines  as 
well  as  strong  vortex  regions.  The  results  were  correlated  by  Hsia  [6] 
and  Hsia  et  al.  [7]  on  the  basis  of  second  order  blast  wave  theory.  Hsia 
found  that  the  data  could  be  correlated  with  a  blast  wave  from  a  line 
charge  parallel  to  the  nozzle  centerline,  having  its  origin  at  the  loca¬ 
tion  of  the  apex  of  the  bow  shock.  A  simpler  correlation  of  shock  shapes 
(and  separation  zones)  was  proposed  by  Wilson  and  Comparin  [8]  on  the 
basis  of  tests  run  for  NASA  by  Vickers  [9].  In  the  Vickers  tests,  shocks 
were  located  by  means  of  interpolation  from  rather  sparsely  populated 
pressure  tap  data  and  are  not  believed  to  be  very  accurate. 


2.  Statement  of  tiie  Problem 

Secondairy  injection  of  a  fluid  into  the  supersonic  flow  in  a 
rocket  nozzle  alters  the  pressure  distribution  on  the  nozzle  wall  con¬ 
siderably  because  of  the  formation  of  a  strong  shock.  Proper  emplacement 
of  injectors  can  thus  be  utilized  to  create  desirable  side  forces.  The 
design  of  such  a  thrust  vectoring  system  requires  accurate  prediction  of 
shock  shapes  and  wall  pressure  distributions  as  functions  of  injection 
location  and  strength.  The  present  study  is  concerned  with  the  first 
problem,  accurate  determination  of  shock  shapes. 

Two  previous  Investigations  of  shock  shapes  in  the  secondary  injec¬ 
tion  in  a  conical  rocket  nozzle  exist;  these  are  the  experiments  by  Charwat 
and  Allegre  [5]  and  the  Vickers  tests  [9] .  Tlae  experiments  of  Chatrwat  and 
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Allegre  were  made  in  5-degree  conical  nozzles  under  a  limited  range  of 
operating  parameters;  those  of  the  Vickers  tests  were  made  in  a  15- 
degree  conical  nozzle.  The  second  order  blast  wave  analysis  of  Hsia 
et  al.  [7]  correlates  the  Charwat  and  Allegre  data,  but  is  at  variance 
with  the  Vickers  tests.  The  present  study  was  made  to  resolve  this  dis¬ 
crepancy  and  to  determine  an  appropriate  approach  to  analysis  of  shock 
shapes  in  15-degree  conical  nozzles. 


3.  Experimental  Approach 

a.  Single  Port  Injection  Studies 

The  shock  experiments  of  this  study  were  performed  as 
part  of  a  larger  set  of  fluidic  secondary  injection  experiments  in  which 
injection  was  accomplished  by  switching  a  bistable  supersonic  valve  into 
opposite  injector  ports.  In  the  present  study,  one  of  the  ports  was 
blocked  with  a  shim  to  eliminate  the  possibility  of  spillover  and  to 
Insure  injection  through  a  single  port. 

Figure  1  schematically  depicts  the  experimental  setup  used  in  this 
study.  A  high  pressure  nitrogen  supply  was  connected  to  a  pressure  regu¬ 
lator  through  a  hand  valve.  The  flow  was  remotely  regulated  by  means  of 
a  fluidic  resistive-capacitive  circuit  which  varied  the  flow  slowly  from 
0  psla  to  a  maximum  of  approximately  1600  psia.  The  supply  flow  feeds 
the  stagnation  chamber  for  the  primary  nozzle  and  secondary  injection 
circuit.  During  the  shock  experiments,  the  solenoid  valve  directed  the 
control  flow  into  the  supersonic  secondary  injection  valve  so  as  to 
direct  its  output  into  the  unblocked  injector  port. 

Two  steel  primary  nozzles  were  used  in  this  stuoy;  dimensions  and 
designations  of  these  nozzles  are  given  in  Figure  2.  The  nozzles  were 
identical  except  for  the  location  of  their  injection  ports.  The  primary 
nozzles  were  bolted  onto  the  cylindrical  stagnation  chamber  and  sealed 
with  a  large  0-ring. 

Twelve  experiments  were  conducted  on  single-port  secondary  injection, 
six  with  each  of  the  two  nozzles  described  in  Figure  2.  Secondary- 
to-prlmary  flow  rate  ratios  were  varied  by  varying  the  nozzle  throat 
diameter  of  the  secondary  injection  valves  (SIVN)  and  the  control  orifice 
diameter.  Test  conditions  are  summarized  in  Table  I.  The  secondary 
throat  diameters  d*  shown  in  Table  1  are  equivalent  diameters  which 

account  for  the  control  flow  through  the  solenoid  into  the  secondary 
injection  flow.  The  runs  listed  in  Table  I  do  not  include  those  experi¬ 
ments  which  yielded  no  usable  shock  data  (runs  1,  2,  7,  and  9). 

The  method  used  to  plot  the  shock  traces  on  the  nozzle  wall  was  that 
first  used  by  Amick  and  Hays  [2].  Before  each  run,  the  Inside  of  the 
nozzle  was  coated  with  a  tacky  mixture  of  china  clay  particles  in  methyl 
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Table  I.  Test  Conditions,  Single  Injection  Port  Studies 


Run 

L 

(in.) 

d* 

s 

(in.) 

n 

3 

0.0976 

0.0450 

4 

mgsM 

0.1496 

0.1058 

5 

0.651 

0.1533 

0.1111 

6 

0.651 

0.1845 

-  0.1609 

8 

0.859 

0.0736 

0.0256 

10 

0.859 

0.0976 

0.0450 

11 

0.859 

0.1030 

0.0501 

12 

0.859 

0.1536 

0.1115 

(PRIMARY  THROAT 
DIAMETER)  d^- 0.46  in 


EQUIVALENT  SECONDARY  /  \ 

THROAT  DIAMETER  '  ' 


salicylate  (Oil  of  Wintergreen)  and  a  thin  layer  of  dry  china  clay  was 
blown  on  top  of  this  mixture.  Along  the  line  where  the  shock  inter¬ 
sected  the  wall,  an  area  was  scrubbed  clean  of  the  china  clay  during 
each  run.  Upon  drying,  this  line  showed  up  as  a  bare  line  on  a  back¬ 
ground  of  white  china  clay.  The  trace  was  plotted  by  measuring  the  dis¬ 
tance  from  the  nozzle  exit  to  the  trace  at  5-degree  intervals  of  azi¬ 
muthal  angle.  To  facilitate  the  data  taking,  a  special  fixture  was 
designed  to  fit  onto  the  exit  of  the  nozzle.  This  fixture  had  a  rotating 
ring  to  measure  and  set  the  azimuthal  angle  to  a  probe  which  could  be 


moved  along  the  wall;  distances  to  the  shock  trace  from  the  nozzle  exit 
were  measured  with  a  machinist's  scale  mounted  on  the  probe.  The  least 
count  on  the  scale  was  0.01  inch  and  the  least  count  on  the  azimuth  ring 
was  5  degrees. 

b .  Adjacent  Port  Injection  Studies 

Two  additional  experiments  were  run  to  determine  the 
shock  shapes  existing  in  a  nozzle  with  injection  from  two  adjacent  ports 
(90-degree  azimuthal  angle  difference).  The  experimental  setup  was 
essentially  the  same  as  that  described  previously  except  that  the  injection 
ports  were  supplied  directly  by  nozzles  instead  of  supersonic  bistable 
valves.  Test  conditions  are  listed  in  Table  II.  The  shock  trace  data 
were  collected  as  described  previously. 


Table  II.  Test  Condicions,  Adjacent  Port  Injection 


d- 

I, 

d* 

P 

(in.) 

d 

Q 

(in 

.) 

Run 

(in. ) 

(in.) 

Bottom  Port 

Left  Port 

13 

0.651 

0.654 

1.363 

0.1442 

0.1448 

14 

0.651 

0.654 

1.363 

0.1442 

0.089 

4.  Analysis  of  Shock  Shapes 

The  formation  of  a  shock  in  a  supersonic  stream  caused  by 
injection  of  a  secondary  fluid  is  analagous  to  the  formation  of  a  bow 
shock  around  an  obstacle.  The  injected  fluid  must  be  turned  downstream 
and  accelerated  through  momentum  interchange  with  the  primary  stream. 
Roughly,  the  penetration  height  of  the  secondary  stream  can  be  related 
to  a  dimension  of  some  equivalent  blunt  body  shape;  the  bow  shock 
around  this  body  may  then  be  studied  and  related  to  the  shock  in  the 
secondary  injection.  This  approach  has  been  used  with  some  success  by 
Hsia,  et  al.  [6|.  The  development  presented  herein  parallels,  with 
slight  modifications ,  the  method  used  by  Narasaki  flO]  in  his  study  of 
the  shocks  around  solid  obstacles  placed  in  a  supersonic  rocket  nozzle. 

Figure  3  depicts  the  shock  shape  in  a  diametral  plane  of  a  nozzle 
with  secondary  injection.  The  X-axis  is  taken  parallel  to  the  nozzle 
centerline  with  its  origin  at  the  intorsectiort  point  of  the  shock  and  the 
nozzle  wall  (upstream  of  the  injection  port).  It  is  assumed  that  the 
shock  shape  is  the  result  of  a  superposition  of  the  primary  nozzle  flow 
and  a  shock  wave  emanating  from  a  line  charge  along  the  X-axis,  the 
strength  of  which  can  bo  related  to  the  size  of  the  equivalent  blunt  body. 
Sakurai  fill  derived  the  following  relation  between  the  radius  R  and  the 
speed  U  of  a  shock  front. 


4 


/  •v2/r 

r  / 

\2 

v4 

(b)(- 

r)  + 

0  "t(v 

- 1 - * 

+ 

whe.-e  and  ^re  constants  and  a  =  0,  1,  and  2  correspond  to 

plane,  cylindrical,  and  spherical  shocks,  respectively.  The  sonic 
velocity  of  the  undisturbed  fluid  outside  the  shock  front  is  represented 
by  c.  The  terra  in  Equation  (1)  is  given  by 


where  E  is  the  energy  release  per  unit  area  of  blast  surface  and  p  is 
the  pressure  of  the  undisturbed  fluid. 

Equation  (1)  applies  to  a  shock  propagating  into  a  quiescent  fluid, 
As  indicated  in  Figure  3,  the  primary  flow  in  a  nozzle  sets  up  a 
standing  shock  with  slope 


s  U  U 
dZ  "  U  ~  cM 


where  is  the  radius  of  the  standing  shock  and  M  is  the  Mach  number  of 

the  primary  (undisturbed)  flow.  Choosing  a  =  1  (cylindrical  shock  front) 
and  using  only  the  second  order  term  in  (c/U),  Equation  (1)  may  be 
written  as 


Rearranging  and  separating  variables. 
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Equation  (4)  may  be  integrated  from  the  shock  apex  (2=0)  where  =  0 
to  a  location  Z  downstream  of  the  shock  apex  to  yield 


r - r~:r  i 

/.w  M  .1  -J-  /■  -ii- 


where  =  denotes  centerline  distance  measured  downstream  from  the  nozzle 
throat  and  corresponds  to  the  location  of  the  shock  apex.  Defining 

f(Z)  by  the  relation 


J 


where  M  denotes  the  free  stream  Mach  number  at  the  shock  apex, 
o 

Equation  (5)  can  solved  for  R^;  thus. 
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Following  Narasaki  (10),’  the  energy  release  per  unit  length  per 
unit  time,  E  is  taken  to  be  exactly  twice  the  drag  on  the  equivalent 
blunt  body.  Thus,  E  =  D/n  and  Equation  (2)  becomes 


R  =  /i.  = 
o  /  p_  / 


where  D  is  the  drag  on  the  equivalent  blunt  body  of  frontal  area  A, 


U  O’ 


6 


In  Equation  (9),  represents  a  drag  coefficient  and  the  subscrint  (  ) 

refers  to  conditions  just  ahead  of  the  shock.  Assuming  the  equivalent 
body  to  be  a  quarter  snhere  with  a  half-cylindrical  afterbody,  Equation 
(9)  may  be  combined  with  Equation  (7)  to  yield,  with  some  manipulation, 


1 

2 

^  ,  (10) 


where  r  represents  the  radius  of  the  quarter  sphere.  Equation  (10) 
differs  from  a  similar  equation  quoted  by  Hsia  et  al  [7]  by  the  factor 
f(Z)  which  accounts  for  the  variation  of  Mach  number  with  Z  and  by  the 
constant  which  has  been  added  here  as  a  constant  to  be  determined 

from  experiments. 


The  drag  coefficient  may  be  estimated  by  considering  the  pressure 
distribution  on  the  nose  of  the  equivalent  blunt  body  (c.f.,  Figure  4). 
As  suggested  by  Zukoski  and  Spaid  [4)  one  may  utilize  Newtonian  flow 
calculations,  as  modified  by  Lees  [12],  to  get  an  expression  for  the 
pressure  coefficient  on  the  surface  of  the  body 


C 

P 


C* 

P 


sin^ 
sin^  B* 


(11) 


where  B  is  the  angle  between  the  local  tangent  to  the  surface  and  the 
undisturbed  primary  flow  and  C  is  defined  by  the  relation 


C 

P 


(12) 


where  the  subscript  (  )  refers  to  conditions  upstream  of  the  bow  shock. 

C*  denotes  the  pressure  coefficient  at  the  point  B  =  B*.  At  the  nose  of 
P 

the  body  B*  =;  rt/2,  and  Equation  (11)  may  thus  be  written 

C  -  C*  sin^  3  .  (13) 

P  P 


7 


Mi 
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The  drag  force  can  be  found  by  integration  of  the  pressure  coefficient 
over  the  surface  of  the  quarter  sphere.  Thus, 


C*'  (X' 


// 


(p- 


p  \  cos  ^  dA 
/  s  s 


(14: 


where  t  denotes  the  angle  between  the  normal  to  the  area  element  and 

the  direction  of  the  free  stream;  A  is  the  projected  area  of  the 
2 

quarter  sphere,  tr  /2.  If  an  area  element  dA^  is  chosen  as  indicated 
in  Figure  5, 

2 

dAg  =  :tr  cos  3  d3 

and 


cos  <>2  =  sin  3 

Thus,  Equation  (14)  may  be  Integrated 


D  = 


Cd 


f  I 

C  CC  I  ^ 

J  2 


Cp  p  U  sin  3  cos  3  ^3 


or 


(15 


The  pressure  coefficient  at  the  nose  of  the  body  may  be  evaluated  as 


C* 

P 


(16 


where  p^  Is  the  stagnation  pressure  behind  a  normal  shock,  easily  calcu' 
lated  or  read  from  tabulated  data  [ 13] . 


The  equivalent  blunt  body  radius  r  is  found  from  an  analysis  of  the 
momentum  interchange  between  primary  and  secondary  fluid  streams. 

Taking  the  case  of  secondary  injection  normal  to  the  nozzle  centerline, 
it  is  assumed  that  the  drag  force  D  is  that  force  required  to  accelerate 
the  injected  fluid  to  the  primary  stream  velocity  as  it  turns  downstream. 
Thus, 


where  m,  is  the  Injectant  flowrate  and  the  subscript  (  )  .  refers  to 

1  ool 

undisturbed  flow  conditions  at  the  injection  station  location.  Equation 
(17)  may  be  manipulated  to  yield 


Upon  substitution  of  Equation  (17)  into  Equation  (10)  the  following 

relation  for  R  results 
s 


Av  A  ,  m. 
«il  i 

itj  ih 

,  0  P  ; 


/  «J  m 


[Zf(z)]‘ 


Equation  (19)  Illustrates  the  curious  fact  that  R  does  not  depend  on 

s 

the  drag  coefficient  C^.  Prestimably,  one  could  use  any  convenient  value; 

this  has  been  pointed  out  previously  by  Hsia  [6] .  However,  selection  of 
does  affect  the  location  of  the  shock  apex  and  thus  its  value  does 

Influence  the  shape  of  the  shock  trace  on  the  wall.  It  can  also  be  shown 
that  Equation  (19)  is  indenendent  of  the  shane  assumed  for  the  equivalent 
blunt  body. 

The  shock  standoff  distance,  from  the  nose  of  the  equivalent  blunt 
body,  is  most  easily  found  by  means  of  the  simple  correlation  suggested 
by  Ambrosio  and  Wortman  [14] 

3.24Al^ 

A  =  0 . 143  r  e  ° 
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Equation  (20)  is  based  on  data  taken  on  spheres  in  supersonic  air  flows 
and,  in  spite  of_its  simplicity,  it  is  relatively  accurate.  Inasmuch 
as  r  depends  on  Cj^,  Equation  (20)  shows  that  A  also  depends  on  C^; 

thus,  one  must  have  a  relatively  accurate  for  an  accurate  estimate 
of  shock  apex  location. 

When  using  Equation  (16),  one  does  not  know,  a  priori,  the  proper 
value  of  M^,  the  Mach  number  at  the  apex  of  the  shock,  :he  location  of 

which  depends  on  C^.  Thus,  proper  calculation  requires  an  iterative 
scheme . 

To  relate  the  shock  shapes  in  a  diametral  plane  to  obseirved  shock 
traces  on  the  nozzle  wall  we  assume  symmetry  of  the  shock  about  the 
Z-axis  (c.f.  Figure  6).  As  Illustrated  in  Fl^jUre  7,  the  location  of 
the  intersection  of  the  shock  with  the  nuzzle  wall  at  each  axial  loca¬ 
tion  is  simply  the  intersection  point  of  two  circles;  one  with  center 
on  the  nozzle  centerline  and  radius  equal  to  the  nozzle  radius  at  that 
location,  the  other  centered  on  the  Z-axis  with  radius  equal  to  the 
shock  radius  at  that  location.  Thus,  at  the  intersection  point,  the 
law  of  cosines  yields 


(21) 


where  a  represents  the  nozzle  radius,  is  the  nozzle  radius  at  the 

shock  apex,  and  is  the  azimuthal  angle  to  the  intersection  point 
measured  relative  to  a  diametral  plane  through  the  centerline  of  the 
injection  port. 

The  preceding  development  completes  the  second  order  blast  wave 
analysis  of  shock  shapes  in  a  nozzle  with  secondary  injection.  Under 
a  given  set  of  conditions,  the  following  procedure  could  be  utilized 
to  apply  the  theory: 

a)  Using  the  primary  flow  Mach  number  at  the  centerline  of  the 
injection  port,  calculate  the  drag  coefficient  C^,  Equation  (16). 

b)  Calculate  the  equivalent  blunt  body  radius  r  from  Equation  (18) 
and  the  shock  apex  location  A  from  Equation  (20). 

c)  Using  the  Mach  mimber  M^  evaluated  at  the  shock  apex  found  in 
.tep  b) ,  Iterate  steps  a)  and  b)  to  the  desired  degree  of  accuracy. 

d)  Apply  Equation  (19)  to  calculate  as  a  function  of  Z. 

e)  Calculate  the  location  of  the  trace  of  the  shock  on  the  nozzle 
wall  from  Equation  (21). 
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Inputs  to  the  procedure  previously  detailed  Include:  primary  and 
secondary  flow  rates,  nozzle  geometry,  and  the  constants  J  from 
Sakural's  solution  (11].  For  y  =  1.4,  a  =  1,  ° 


=  -1.989 


Jq  =  0.877 


In  evaluating  as  a  function  of  Z  (step  d  cf  the  preceding  pro¬ 
cedure)  ,  from  Equation  (19) ,  It  Is  necessary  to  evaluate  the  function 
f(Z).  Equation  (6)  defines  this  function  as  an  Integral  involving  the 
free  stream  Mach  number.  The  variation  of  free  stream  Mach  number  with 
distance  is  plotted  in  Figure  8  for  each  of  the  primary  nozzles  used 
In  this  study.  It  was  found  that,  with  the  exception  of  a  short  region 
near  the  throat,  a  good  approximation  to  the  Mach  number  variation  Is 
given  by  an  expression  of  the  form 


The  specific  expressions  used  In  this  study  are  compared  to  the  actual 
Mach  number  variations  in  Figure  8. 


5.  Results  and  Oiscunion 

a.  Shock  Shapes,  Single  Port  Injection 


Data  on  shock  patterns  were  gathered  (as  described  In 
Section  3)  in  a  series  of  twelve  experiments  conducted  on  two  15-degree 
conical  nozzles.  Both  nozzles  had  Identical  throat  and  exit  diameters 
of  0.460  and  1.155  Inches,  respectively;  they  differed  only  in  the 
location  of  their  injection  ports  (see  Figure  2  and  Table  I  for  descrip¬ 
tion  and  dimensions) .  The  experiments  covered  a  range  of  secondary- 
to-primary  flow  rate  ratios  from  2.5  to  16  percent.  Nitrogen  was  used 
as  the  working  fluid  in  all  experiments. 

Figures  9  through  12  depict  the  shock  shapes  for  single  port  injec¬ 
tion  into  the  first  nozzle  (designation:  SIV003)  which  had  its  injection 
port  approximately  midway  between  the  throat  and  exit.  The  data  for 
the  second  nozzle  (SIV004) ,  with  injection  port  approximately  one-third 
the  nozzle  length  downstream  of  the  throat,  is  illustrated  in  Figures 
13  through  16.  In  both  cases,  the  data  are  presented  in  the  order  of 
increasing  injected  flow  rate.  All  of  the  cuirves  shown  are  plotted  on 
grids  representing  the  unwrapped  nozzle  walls. 
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Several  features  were  connon  to  all  the  tests.  First,  considerable 
data  loss  was  experienced.  Thus,  e.g.,  in  Run  3  (Figure  9)  only  a  short 
section  of  the  china  clay  streak  pattern  was  observed  on  one  side  of  the 
injection  port;  presumably,  the  remainder  of  the  china  clay  particles 
were  blown  off  before  the  completion  of  the  run.  Coating  of  the  nozzle 
with  the  china  clay-methyl  salicylate  mixture  before  each  run  remained 
somewhat  of  an  art  throughout  the  experiments  and  data  loss  was  experi¬ 
enced  to  some  extent  in  all  runs.  Most  of  the  wall  traces  consisted  of 
bare  areas  on  white  backgrounds  of  the  coating  of  china  clay  particles. 
Where  discernible,  the  downstream  edge  of  th.e  bare  area  was  more  dis¬ 
tinct  than  the  upstream  edge.  In  one  case,  the  upstream  edge  was  too 
faint  to  be  accurately  located  (Run  6,  Figure-  12);  in  another  (Run  8, 
Figure  13)  the  rear  edge  of  the  trace  could  not  be  located,  apparently 
because  the  china  clay  particles  were  blown  away  during  the  run.  In 
those  runs  where  both  edges  of  the  trace  could  be  accurately  located 
(e.g..  Run  5,  Figure  11)  the  upstream  and  downstream  edges  tend  to  merge 
as  the  shock  wraps  around  the  wall.  If  the  upstream  edge  is  interpreted 
as  the  separation  zone  ahead  of  the  shock  at  the  rear  edge  of  the  trace, 
the  merging  of  the  two  is  anticipated  in  agreement  with  the  results  of 
Charwat  and  Allegre  [5]. 

Comparisons  of  measured  shock  traces  with  those  predicted  by  second 
order  blast  wave  theory  (presented  in  Section  4)  are  given  in  Figures  9 
through  16.  In  general,  overall  agreement  is  fair  if  the  rear  edge  of 
the  measured  trace  is  taken  to  be  the  bow  shock  and  the  constant  Cj^  is 

taken  as  1.33  (~  At  the  lower  flow  rate  ratios,  application  of 

the  theory  is  complicated  by  the  size  of  the  injection  port  (0.287-inch 
diameter  in  each  case).  The  Mach  number  of  the  secondary  fluid  at  the 
injection  port  is  unknown;  if  it  were  assumed  to  be  unity,  the  injectant 
stream  would  never  completely  fill  the  injection  port.  The  procedure 
adopted  was  to  take  the  injectant  Mach  number  as  unity  and  to  consider 
the  injectant  strea  to  be  exiting  from  the  rear  of  the  injection  port. 

The  shock  standoff  distance  was  then  calculated  relative  to  the  center- 
line  of  that  stream.  In  Runs  3,  8,  10,  and  11  this  procedure  resulted 
in  calculated  shock  apexes  inside  the  Injection  port.  For  those  cases 
the  shock  apex  was  arbitrarily  taken  at  the  upstream  edge  of  the  injection 
port,  resulting  in  good  agreement  between  measured  and  predicted  shock 
shapes  for  Runs  10  and  11;  Run  3  was  not  accurately  predicted  with  this 
technique  and  Run  8  yielded  no  rear  edge  with  which  to  compare. 

Since  Cj >  1  for  quantitative  agreement  between  theory  and  experiment, 

it  can  be  concluded  that,  near  the  wall,  the  shock  spreads  more  rapidly 
than  theory  predicts.  It  is  believed  that  this  behavior  is  peculiar  to 
the  present  geometry  and  is  not  to  be  extrapolated  to  other  geometries . 

Shock  apex  location  was  accurately  predicted  for  Runs  4  and  5  by 
using  the  Mach  number  of  the  free  stream  at  the  center  of  the  injection 
port  area,  assuming  sonic  Injection,  and  performing  the  iteration 
described  in  Section  4.  However,  at  the  higher  flow  rate  ratios  (Runs 
6  and  12)  the  location  of  the  shock  apex  is  drastically  over  predicted  by 
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theory  and  no  rational  method  of  accurately  predicting  shock  apex  location 
was  found;  the  theoretical  curves  were  forced  to  pass  through  the 
measured  apexes  in  those  cases. 


Insufficient  simultaneous  data  on  separation  and  bow  shocks  exist 
to  get  a  correlation  between  the  two.  Run  5  indicates,  however,  that 
the  two  are  separated  at  the  shock  apex  by  approximately  1.5  r.  The 
Wilson  and  Comparin  18]  estimate  of  this  distance  is  1.75  r;  their 
correlation  of  the  separation  line  is 


(f)  =  "=(?) 


where  x  denotes  distance  along  the  wall  from  the  apex  of  the  separation 
line.  This  relation  predicts  a  separation  line  which  diverges  from  the 
bow  shock,  whereas  measurements  indicate  a  merging  of  the  two.  No 
attempt  was  made  to  Improve  this  correlation. 


The  modification  of  second  order  blast  wave  theory  to  account  for 
a  varying  free  stream  Mach  number  results  in  a  slight  improvement  as 
illustrated  in  Figure  11  for  Run  5.  As  indicated,  if  no  correction 
is  made  for  variation  of  free  stream  Moch  number  (i.e.,  f  «  1),  the 
theory  drastically  overpredicts  the  shock  spread.  For  the  conditions 
of  Run  5,  in  fact,  shock  impingement  on  the  opposite  nozzle  wall  is 
predicted  but  not  measured. 


Wilson  and  Comparin  [8]  proposed  the  following  correlation  of  shock 
shape  .s ; 


(?)  “  (?) 


0.535 


where  (a<i)  represents  the  dlstanc<j  measured  along  the  wall  from  a  diame¬ 
tral  plane  passing  through  the  centerline  of  the  injection  port  and  x* 
is  the  distance  along  the  wall,  in  that  diametral  plane,  measured  from 
the  shock  apex.  As  Illustrated  in  Figure  11,  Equation  (25)  drastically 
underpredicts  the  spread  of  the  shock.  It  was  found  that  at  the  lower 
flow  rates  (less  than  10  percent  of  the  primary  flow)  a  fair  correlation  of 
shock  shapes  could  be  obtained  if  the  constant  1.12  in  Equation  (25)  is  changed 
to  2.13.  A  number  of  these  curves  are  .shewn  in  Figures  9  through  12. 


Shock  Shapes,  Adjacent  Port  Injection 


Two  experiments  were  conducted  with  Injection  from  two 
adjacent  ports  (90-degrees  azimuthal  angle  separation) .  These  experi¬ 
ments  were  run  on  a  15-degree  conical  nozzle  with  a  0.654-inch  diameter 
throat  and  a  1.363-lnch  diameter  exit  (designated  STV006) .  The  first-run 
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was  made  with  vxearly  identical  injection  nozzles  (designated  f?IVN002-l, 
SIVNOOl-2) ;  the  second  run  was  made  with  injection  nozzles  of  signi' 
ficantly  different  throat  diameters  (SIVN006,  SIVN002-2).  A  complete 
listing  of  all  geometry  is  shown  in  Table  II. 


The  results  of  the  first  run,  with  nearly  identical  nozzles  (i.e., 
a  syi.unetrical  shock  pattern) ,  is  shown  in  Figure  17  together  with  the 
shock  shapes  predicted  with  second  order  blast  wave  theoiry.  Agreement 
between  the  measurements  and  predictions  is  good.  Downstream  of  the 
intersection  point  (at  45  degrees)  is  a  slip  line  extending  to  the  exit 
of  the  nozzle.  Shock  reflections  would  be  anticipated  for  such  an 
intersection,  but  these  were  not  detected  in  the  present  series  of 
measurements. 


Measured  and  predicted  shock  shapes  for  the  second  run,  Tjith  an 
asymmetrical  pattern,  are  given  in  Figure  18.  Again,  the  agreement  is 
good.  In  this  case  the  shock  intersection  occurs  at  52  degrees  left 
of  the  injector  port,  at  the  bottom  of  the  nozzle,  yet  the  slip  line 
runs  straight  hack  to  the  exit.  In  both  tests  the  slip  line  downstream 
of  the.  intersection  point  was  exceptionally  clear. 

In  Figures  17  and  18,  the  predicted  shock  shapes  based  on  second 
order  blast  wave  theory,  were  forced  to  pass  through  the  measured  apexes 
of  the  shock  traces.  The  shock  apexes  were  again  slightly  overpredicted 
(0 . 09  inch) . 

An  important  bonus  was  derived  from  the  shock  intersection  tests. 
Side  force  measurements  were  made  with  each  port  off,  each  port  on, 
and  both  ports  on.  Tliese  measurements  indicate  zero  coupling  of  side 
forces,  even  though  the  bow  shocks  intersected.  Tliese  results  agree 
with  the  force  measurements  of  the  Vickers  tests  [9]  . 
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All  of  the  data  previously  presented  were  found  to  be  reproducible 
to  within  0.01  inch,  the  least  count  on  the  scale  used  to  position  the 
probe . 


6.  Conclusions 

The  results  of  this  study  lead  to  the  following  conclusions 
concerning  shocks  in  secondary  injection  into  15-degree  conical  nozzles: 

a)  The  shape  of  the  shock  trace  on  the  nozzle  wall  is  well  repre¬ 
sented  by  a  second  order  blast  wave  analysis,  which  accounts  for  varia¬ 
tion  of  free  stream  Mach  number  along  the  nozzle  centerline  except  for 
those  cases  where  the  injection  ports  do  net  flow  full.  Good  quantita¬ 
tive  agreement  between  theory  and  experiment  was  obtained  by  proper 
choice  of  the  empirical  constant  (r  .^jT  for  the  present  geometry) . 

b)  A  separation  shock  exists  a  short  distance  upstream  of 
the  shock  trace  on  the  wall  (within  1.5  r).  This  separation  shock  merges 
with  the  bow  shock  as  it  wraps  around  the  wall  of  the  nozzle. 
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Measurements  indicate  a  slip  line  extending  straight  downstream  from 
this  point  to  the  exit,  even  with  unequal  strength  shocks. 

d)  The  side  force  produced  by  secondary  injection  from  one 
port  is  completely  unaffected  by  injection  from  an  adjacent  (90  degrees) 
port  with  or  without  bow  shock  intersection. 


Figure  8.  Mach  Number  Versus  Axial  Distance  from  Throat 


'  '  if' 


Ferrari,  C.,  Interference  Between  a  Jet  Issuing  Laterally  from  a 
Body  and  the  Enveloping  Supersonic  Stream.  The  John  Hopkins  Univer¬ 
sity,  Applied  Physics  Laboratory,  April  1959,  Bumblebee  Report  No. 
286. 

Amick,  J.  L.  and  Hays,  P.  B.,  Interaction  Effects  of  Side  Jets 
Issuing  from  Flat  Plates  and  Cylinders  Aligned  with  a  Supersonic 
Stream,  Wright  Air  Development  Division,  June  1960,  Report  No. 

TR  60-329. 

Broadwell,  James  E.,  "Analysis  of  the  Fluid  Mechanics  of  Secondairy 
Injection  for  Thrust  Vector  Control,"  AlAA  Journal.  1^,  No.  5, 

May  1963,  pp.  1067-1075.  (See  also  Broadwell,  James  E.,  "Corre¬ 
lation  of  Rocket  Nozzle  Gas  Injection  Data,"  AIAA  Journal,  1,  No.  8, 
August  1963,  pp.  1911-1913.) 


/ 


4.  Zukoskl,  Edward  E.  and  Spaid,  Frank  W. ,  "Secondary  Injection  of 
Gases  into  a  Supersonic  Flow,"  AIAA  Journal,  2,  No.  10,  October 
1964,  pp.  1689-1696. 

5.  Charwat,  A.  F.  and  Allegre,  J.,  "Interaction  of  a  Supersonic  Stream 
and  a  Transverse  Supersonic  Jet,"  AIAA  Journal,  2,  No.  11,  November 
1964,  pp.  1965-1972. 

6.  Hsia,  Heniry  T-S.,  "Equivalence  of  Secondary  Injection  to  a  Blunt 
Body  in  Supersonic  Flow,"  AIAA  Journal,  4,  No.  10,  November  1966, 
pp.  1832-1834. 

7.  Hsia,  H.  T-S.,  Seifert,  H.  S.,  and  Karamcheti,  K.,  "Shocks 
Induced  by  Secondary  Fluid  Injection,"  Journal  Spacecraft.  2,  No.  1, 
January- February  1965,  pp.  67-72. 

8.  Wilson,  W.  G.  and  Comparin,  R.  A.,  "Analysis  of  the  Flow-Disturbance 
and  Side  Force  Due  to  Gaseous  Secondary  Injection  into  a  Rocket 
Nozzle,"  Journal  Spacecraft.  2>  No.  5,  May  1970,  pp.  539-543. 

9.  Vickers,  Inc.,  "Proportional  Solid  Propellant  Secondary  Injection 
Thrust  Vector  Control,"  November  1966,  NASA  CR-637. 

10.  Narasaki,  T.,  "The  Side  Force  Induced  by  an  Obstacle  in  the  Nozzle 
of  a  Rocket  Engine,"  Paper  presented  to  the  XXI®^  international 
Astronautical  Congress,  Korstanz,  Federal  Republic  of  Germany, 

4"10  October  1970. 

11.  Sakurai,  A.,  "On  the  Propagation  and  Structure  of  a  Blast  Wave," 
Journal  Phys.  Soc.  Japan.  Part  I;  8,  No.  5,  September-October  1953, 
pp.  662-669;  Part  II:  9,  No.  2,  March-April  1954,  pp.  256-266. 


31 


d 


12.  Lees,  L.,  "Hypersonic  Flow,"  IAS  Paper  554,  1955. 


13.  Ames  Research  Staff.  Equations,  Tables,  and  Charts  for  Compressible 
Flow,  1953,  NACA  1135. 


14.  Ambrosio,  A.  and  Wortman,  A.,  "Stagnation  Point  Shock  Detachment 
Distance  for  Flow  Around  Spheres  and  Cylinders,"  ARS  Journal,  32 , 
February  1962,  p.  281. 


32 


